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ABSTRACT 

We investigate the possible contribution from the emission of accretion flows 
around supermassive black holes in early type galaxies to current measurements of 
the Cosmic Microwave Background (CMB) at radio frequencies. We consider a range 
of luminosities suggested by targeted radio observations and accretion models and 
compute the residual contribution of these sources to the spectrum and bispectrum of 
the observed CMB maps. As for high-resolution CMB experiments, we find that the 
unresolved component of these sources could make up to ~ 40 — 50% of the observed 
CBI and BIMA power spectrum at I > 2000. As a consequence, the inferred erf z 
value could be biased high by up to 6 — 7%. As for all sky experiments, we find that 
the contribution of accretion-flow sources to the WMAP bispectrum is at the 2-3 per 
cent level at most. At the flux limit that Planck will achieve, however, these sources 
may contribute up to 15 per cent of the bispectrum in the 60-100 GHz frequency 
range. Moreover, Planck should detect hundreds of these sources in the 30-300 GHz 
frequency window. These detections, possibly coupled with galaxy type confirmation 
from optical surveys, will allow number counts to put tighter constraints on early-type 
galaxies radio luminosity and accretion flows properties. These sources may also con- 
tribute up to the 30 per cent level to the residual radio sources power spectrum in 
future high-resolution SZ surveys (like ACT or APEX) reaching mJy flux limits. 

Key words: accretion - cosmology: cosmic microwave background - galaxies: ellip- 
tical - galaxies: nuclei - infrared:galaxies 



1 INTRODUCTION 

Cosmic microwave background (CMB) experiments have 
been able to measure cosmological parameters to an un- 
precedented level of accuracy (Spergel et al. 2003). For this 
to be possible, however, the contribution to the primary sig- 
nals must be disentangled from other astrophysical emissions 
at the observed frequencies. Therefore, foreground identifi- 
cation and removal is fundamental. 

One of the major sources of contamination at small 
scales is constituted by point sources, and several studies 
have been carried out to estimate their contribution to CMB 
anisotropy experiments (Toffolatti et al. 1998; White & Ma- 
jumdar 2004; Argiieso, Gonzalez-Nuevo & Toffolatti 2003). 
Generally, source counts are determined in the radio band, 
through deep VLA surveys down to /J,3y levels at 1.41, 4.86 
and 8.44 GHz. These counts are then extrapolated to the 
higher range of frequencies relevant for CMB experiments 
(Toffolatti et al. 1998; Toffolatti & et al. 1999; De Zotti et al. 



2000; Gawiser & Smoot 1997). This extrapolation provides 
a reasonable estimate of the contribution from the "steep" 
and "flat" spectrum sources (with F v cx v~ a , and a ^ 0, 
such as compact radio galaxies and radio loud quasars), but 
it strongly under-represents an important contribution from 
a class of sources with inverted spectra (a < 0; e.g. De Zotti 
et al. 2000). 

Inverted-spectrum sources, such as GHz peaked sources 
(GPS; O'Dea (1998); Guerra, Haarsma & Partridge (1998)), 
bright (flux ~ 1 — 10 Jy) and rare, are generally associated 
with bright active galaxies or quasars at high redshifts. Their 
contribution to the CMB experiments was studied by De 
Zotti et al. (2000). These sources, being bright and rare, are 
easily identified in the CMB maps, and removed. 

There is however another class of inverted-spectrum 
sources, much fainter (flux ~ 1 mJy) but much more com- 
mon, associated with emission from the nuclei of nearby 
galaxies. Radio continuum surveys (at v < 8 GHz) of ellip- 
tical and SO galaxies have shown that the sources in radio- 
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quiet galaxies tend to be extended but with a compact com- 
ponent with relatively flat or slowly rising radio spectra. Re- 
cent VLA studies at high radio frequencies (up to 43 GHz), 
although carried out only on a limited sample of objects, 
have shown that all of the observed compact cores have 
spectra rising up to ~ 20 — 30 GHz. These sources of radio 
emission are believed to result from the process of accretion 
of gas into the supermassive black holes likely ubiquitous in 
the centers of galaxies (Magorrian et al. 1998). 

A popular model to describe the broadband spectral 
energy distribution from the accretion flows around these 
supermassive black holes is the advection-dominated accre- 
tion flow (ADAF) model (Rees et al. 1982; Narayan & Yi 
1994). Within the context of this model, the foreground 
contribution to the CMB experiments, and in particular 
to the Planck Surveyor, was studied by Perna & Di Mat- 
teo (2000). In particular, they estimated the contribution 
of these sources both within the "standard" ADAF model, 
where the accretion rate is independent of the radius R of 
the flow (implying that all the mass is accreted into the 
black hole (BH)), and within the context of the ADAF model 
with winds, where not all the mass is accreted by the BH, 
but some is lost into winds (Blandford & Begelman 1999; 
Narayan & Yi 1995a). In this case the radio emission is sup- 
pressed because winds remove mass from the inner regions 
of the flow, where the synchrotron radiation giving rise to 
the radio emission is produced. Detailed modeling has only 
been possible for a handful of sources so far, as it requires 
multi wavelength, high-resolution data. Observations have 
yielded somewhat mixed evidence: in some cases the emis- 
sion is consistent with ADAFs with winds (i.e. suppressed 
radio emission; Di Matteo et al. 2000), while in others it is 
higher than the standard model with no wind (Di Matteo 
et al. 2001), probably because of extra power output by jets 
associated with the accretion flow. 

While studies of this class of sources of radio emission 
so far have been limited by low number statistics, Perna & 
Di Matteo (2000) showed that CMB experiments with the 
future Planck surveyor mission can provide interesting con- 
straints from a statistical point of view. At the same time, 
these sources can be an important contaminant to CMB 
maps. 

In this paper, we compute the contribution of the emis- 
sion from accretion flows in early-type galaxies to the current 
signal of CBI, BIMA, WMAP. We show that the signal from 
these unresolved sources is able to influence current limits 
on the power spectrum normalization as- Finally, we make 
predictions on the power spectrum from residual sources in 
future high-resolution CMB experiments and on the bispec- 
trum in all-sky experiments like Planck and WMAP. 



2 MICROWAVE EMISSION FROM 

ACCRETION FLOWS IN ELLIPTICAL 
GALAXIES 

As discussed above, there is evidence for the existence of su- 
permassive massive black holes (BHs) at the center of galax- 
ies. Inferred black hole masses appear to be proportional 
to the mass of the bulge component of their host galaxies. 
Therefore, central BH masses are expected to be much larger 
in elliptical than in spiral galaxies. Independently of the de- 



tails of the model for the production of the radio/microwave 
emission, this radiation is expected to be some fraction of 
the accretion energy, which is proportional to the mass of 
the black hole. Therefore the contribution from the nuclei 
of the ellipticals is expected to dominate that from the nu- 
clei of the spirals, which we neglect here. Similarly, among 
the population of ellipticals, the main contribution derives 
from the most massive galaxies. To make a conservative es- 
timate of the emission from these galaxies, we consider only 
the bright end of the distribution, with L Z L*. Studying 
a sample of nearly 9000 ellipticals in the Sloan Digital Sky 
Survey, Bernardi et al. (2003) estimated a comoving num- 
ber density $» = (5.8 ± 0.3) x 10" 3 /i 3 Mpc -3 . The redshift 
evolution of their sample appeared consistent with the law 
$*(z) = 10°- 4Pz $»(0) (P » -2), found by Lin et al. (1999) 
in their sample of galaxies drawn from the Canadian Net- 
work for Observational Cosmology Field Galaxy Redshift 
Survey. 

For a cosmological population of sources with intrinsic 
luminosity function /(L^) and redshift evolution $»(z), the 
differential number counts are given by 



dn(S v ) 
dS v 



dL v 



f(L„)9(z] 



dV{z) 
dz 



dz 
dS v 



(1) 



where dV(z)/dz is the comoving volume. In a flat cosmology 
with a cosmological constant it is given by 
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We assume a cosmological model with fl m =0.3, Qa — 0.7 
and H = 71 km s" 1 Mpc" 1 (Spergel et al. 2003). In Eq.Q, 
z(S v ,L v ) is derived by inverting the relation S u = L v (l + 
z)/A-kD\{z), where Dl(z) is the luminosity distance. 

We will consider two models for the emission from 
accretion flows in early- type galaxies (see figure 0: while 
these are computed on the basis of a specific accretion 
model, they can be considered generally representative of 
the typical range in luminosity and spectral shape that 
characterize these sources. However, we should emphasize 
that, given the small sample of observations so far, and the 
uncertainties in the specific model parameters discussed 
below, our conclusions should be taken as indicative rather 
than general and final. Tighter constraints (which will leave 
less space for parameter variations) will be obtained with 
future, more sensitive CMB experiments together with a 
detailed, individual study of a larger sample of sources. 

a) Model A: 

This is the standard ADAF model (Narayan & Yi 
1994), where the accretion rate is a constant function of the 
radius within the flow. In this model, the radio/microwave 
emission is due to synchrotron emission from the inner 
regions of the accretion flow. The emission at the self- 
absorbed synchrotron peak scales as L v oc v%T, where 
u c oc T 2 B oc T 2 M 1/2 M^R- 5/i (Narayan & Yi 1995a); 
T is the electron temperature, and B the magnetic field 
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strength. The ADAF model depends on a number of 
microphysical parameters (i.e. the flow viscosity, the ratio 
of gas to magnetic pressure, the adiabatic index of the fluid, 
the fraction of turbulent energy which goes into heating 
the electrons) as well as geometrical parameters (i.e. the 
radial extent of the flow, the possible transition from a 
hot state to a cool disk at some radius), as well as on 
the BH mass Mbh and on the accretion rate M. Whereas 
some constraints on these parameters can be derived by 
multiwavelength observations and broadband modeling of a 
given source, they cannot be well constrained a priori and 
they can in principle vary from source to source. Therefore, 
in a statistical study like ours, we can only adopt "typical" 
values for all these parameters. In particular, we use the 
ones adopted by Di Matteo et al. (2000) for a BH of 
mass ~ 1O 9 M0 (as typical of the bright end of the galaxy 
luminosity function that we are considering here) and for 
an accretion rate ~ 0.005 in Eddington units (note that the 
transition to an ADAF is believed (Narayan & Yi 1995b) 
to occur at rates below ~ a 2 MEdd, a being the viscosity 
parameter of the flow). The luminosity predicted by this 
model is shown by the upper line in Figure 1. To allow for 
a spread in BH masses, accretion rates, and parameters of 
the flow around this typical value, we take the luminosity 
function f{L v ) to be a log-gaussian with mean given by the 
curve displayed in the figure, and a standard deviation a. 
We adopt a = 0.25; the predicted number counts are very 
little dependent on the precise value of a. 

b) Model B: 

The luminosity function for this model is calibrated on 
the sample of those sources for which the broadband energy 
distribution was best fit with an ADAF model with winds 
(Di Matteo et al. 2000). At frequencies in the range ~ 30-100 
GHz, the mean luminosity is roughly constant, and on the 
order of L„ ~ 10 28 erg s _1 Hz -1 (see the dashed region in 
figure^ which encompasses the range of those observations). 
Hence, for this model we assume the mean luminosity to be 
constant over the considered frequency range and take the 
luminosity function f(L„) to be a log-gaussian with mean 
given by the above value, and a standard deviation a — 0.25. 

Figure shows the number counts predicted by the 
sources in Model A and B. These sources are typically fainter 
than the radio sources presented in Toffolatti et al. (1998), 
and therefore the number counts are dominated by the local 
population. As a consequence, the slope in figure [5] is very 
close to Euclidean down to low flux levels. Being faint, these 
sources are not likely to be individually detected. However, 
because of the steepness of their number counts, they may 
contribute significantly to the residual signal in CMB exper- 
iments, where their relevance with respect to the Toffolatti 
et al. (1998) type of sources becomes increasingly impor- 
tant as the detection flux limit of the survey is reduced. 
Moreover, because of the radio spectrum of these sources 
is inverted, they are less likely to be found at low frequen- 
cies, where most point sources catalogs are compiled, and 
therefore may be neglected in present data analysis. In the 
following, we will estimate the contribution of these accre- 
tion flow sources to current CMB experiments, and explore 
the consequences that the emission from these (mostly un- 
resolved) point sources can have on the determination of the 
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Figure 1. Synchrotron emission from low-radiative efficiency ac- 
cretion flows. The uppermost curve is the standard ADAF model 
with no outflows, while the shaded region in between the two 
lower curves represent a range of ADAF models with large out- 
flows. The luminosity function in our model A is a log-gaussian 
with mean given by that curve. In the range of frequencies 30-100 
GHz that we are mostly interested in here, the average luminosity 
in the shaded region is ~ 10 28 erg/s/Hz. Our model B is described 
by a log-gaussian with mean luminosity given by this value. 

matter power spectrum normalization as from these exper- 
iments. 



3 POINT SOURCES AND CMB 
EXPERIMENTS 

Point sources are detected in CMB maps as high point-like 
fluctuations above the mean. Their detection depends on in- 
strumental properties such as the beam size and the noise 
level, as well as on the characteristics of the competing sky 
signals. Assuming that all sources above a given flux limit 
Sum are subtracted from the data, residual sources still con- 
tribute to the statistics of the map (Pierpaoli 2003). For 
Poisson-distributed sources, the residual power spectrum 
and bispectrum are constant for all scales and read: 

fS lim 

a 2 {v) = Ci(y) = g(x) 2 / dS dn/dS S 2 (4) 
Jo 

fSlim 

b{u) = g(xf / dS dn/dS S 3 (5) 
Jo 

where g(x) = 2{hc) 2 /(kT) 3 [sinh(x/2)/x 2 ] 2 (with 
(hc) 2 /(kT) 3 = 0.02 M Ksr/Jy) and x = ^(GHz)/56.78. 

The bispectrum is sometimes quoted as the adimen- 
sional quantity b = b/T 3 . If the sources are clustered, the 
power spectrum is increased by: 

Ci{v) clu =w t (I(u)f (6) 
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Figure 2. The number counts of accretion flow sources at differ- 
ent frequencies. Long-dashed (short-dashed) lines correspond to 
model A (B), as described in the text. The solid line are radio and 
infrared sources as reported from Toffolatti ct al. (1998), rescaled 
by a factor 0.8 (Argueso et al. 2003). Fluxes are in Jy, number 
counts in deg~ 2 . 

where the intensity I(v) = g(x) J^ l ' m dS dn/dS S and the 
wi are the Legendre expansion of the angular correlation 
function: 



w(9) = wi Pi (cost 



(7) 



In order to compute the clustering term, we use the recent 
determination of the SDSS red sample angular correlation 
function (Budavari et al. 2003). This estimate should pro- 
vide a conservative limit to the clustering contribution, since 
the red sample contains more galaxies than just ellipticals, 
and ellipticals are more clustered than other sources. 

Finally, given a particular experiment it is useful to 
quote the (non-clustered) power spectrum as the noise per 
pixel caused by the residual point sources as: 



(8) 



where 9 is the pixel dimension, typically assumed as the 
FWHM of the beam. 

In the following, we will estimate the contribution of 
these sources to the current WMAP, CBI and BIMA signal. 

The WMAP spectrum is clearly dominated by the pri- 
mary CMB and, given the high flux cut (~ Uy), the ob- 
served bispectrum is likely to be dominated by other radio 
sources more luminous than elliptical galaxies at these fre- 
quencies. WMAP, however, presents the advantage of being 
an all-sky survey. Number counts of point sources which are 
confirmed to be elliptical galaxies by mean of other (optical 
and infrared) observations may give a hint on the typical 
radio emission of early-type galaxies. 

CBI and BIMA, on the contrary, observed a small area 
of sky with high resolution and sensitivity. These experi- 
ments detected the CMB power spectrum at very high I's, 
where the primary CMB is weak and the point source con- 



tribution may be dominant. Therefore, the contribution of 
the emission from accretion flows in early-type galaxies may 
be significant to the observed CBI and BIMA power spec- 
trum. We will estimate the importance of this contribution 
and its consequences on the derived cosmological parameter 
estimation. 



3.1 WMAP detected sources 

The WMAP experiment (Bennett et al. 2003) detected 
208 point sources with a flux S > 0.75 in the V band. 
Among those sources, 29 have been identified as galaxies 
while 5 have no identification and could therefore be galax- 
ies (Trushkin 2003). Hence at most there is a total of 34 
sources which are candidate for harboring a source of emis- 
sion from an accretion flow. This classification, however, 
does not distinguish elliptical galaxies from spirals. In or- 
der to determine how many of the WMAP sources could 
be elliptical emitting according to the model described here, 
we proceeded as follows. We selected the sources in Trushkin 
(2003) according to the following two criteria: i) we elimi- 
nated sources with measured redshift greater than z ~ 0.01, 
because if they were ellipticals emitting according to our 
models they would have an observed flux lower than 1 Jy; 
ii) we kept only sources with flat or inverted spectra. We 
then visually inspected the remaining sources in the 2MASS 
database, in order to select those which appear to be el- 
liptical galaxies. We only found two sources of this kind. 
A comparison with the SDSS data (which currently covers 
about 10 per cent of the sky) finds only one elliptical galaxy 
among the WMAP sources. A total number of 2-3 ellipticals 
in the WMAP sample of detected sources seems a conserva- 
tive estimate. 

Model A predicts about 4.5 (8) sources above a flux 
Sum = 1(0.75) Jy at 44 GHz, while model B would pre- 
dict 0.15 (0.25) sources for the same flux cut. These predic- 
tions seem to favour model A, model B predictions being 4-6 
a away from observations. Current data therefore appears 
to suggest that the bulk of the elliptical galaxy population 
emits at ~ 40 GHz more than the ~ 10 28 erg s _1 Hz -1 level 
assumed in model B. 

Planck should reduce the source detection threshold of 
at least one order of magnitude with respect to WMAP, pos- 
sibly allowing for the detections of hundreds of these sources. 
Moreover, SDSS at completion should observe about one 
quarter of the sky. By comparing these two data sets it will 
be possible to infer more conclusive statements on accretion 
models purely by means of radio number counts. 



3.2 Current limits from CBI and BIMA 

Recent results on the CMB power spectrum from high reso- 
lution CMB experiments around 30 GHz like CBI and BIMA 
have shown a rise in the high i's which is interpreted as sig- 
nature of the SZ power spectrum from galaxy clusters (Ko- 
matsu & Seljak 2002; Bond 2002). This interpretation allows 
to infer a value of the matter power spectrum as — 1 (Read- 
head 2004). Residual radio point sources, however, are also 
likely to contribute to the observed power spectrum. We ar- 
gue here that accretion flow sources may present a sizable 
contribution to the CBI and BIMA power spectra. 
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In CBI and BIMA individual point sources have been 
subtracted from the data when identified by means of other 
observations of the same sky area at much lower frequencies. 

In the CBI case, point sources above the flux limit of 
3.4 mJy in the NRAO VLA Sky Survey (NVSS) at 1.4 GHz 
have been subtracted from the data. A residual point source 
signal, compatible with the observed point source popula- 
tion above the flux cut, has been modeled in the data anal- 
ysis. The population of sources considered in our work, how- 
ever, is not likely to be well represented by the detected one, 
because our sources are typically very faint at such low fre- 
quencies. As a result, the observed population may be rather 
biased toward flat or falling spectra sources. The observed 
population, in fact, shows a less steep dn/ds (oc 5 1-1 ' 875 ) 
than the one appropriate for the point sources considered 
here ( oc S 1-2 ' 5 ), and a mean spectral index a ~ —0.45, 
while our sources have a — 0.5 — 0.8 in the relevant fre- 
quency range of ~ 1-10 GHz. We extrapolate the NVSS flux 
limit to 30 GHz using a conservative spectral index a — 0.5, 
and obtain for CBI Snm^oGHz = 16 mJy as a nominal flux 
limit for our sources. 

In the case of BIMA, sources are detected with the VLA 
at 4.8 GHz with a flux limit of 150-175 fiJy (Dawson et al. 
2002). If radio sources have a flat or falling spectrum, then 
this procedure ensures the subtraction of sources at 30 GHz 
down to a very low flux limit. For this reason, the BIMA 
data analysis did not consider a possible contribution from 
residual point sources. 

Because accretion flow sources present an inverted spec- 
trum, we do expect a residual signal from them below the 
flux limit at 30 GHz. Considering the conservative case of a 
spectral index of 0.5, such limit is Siim,30GHz = 0.4mJy 

Given the considerations above, we assumed that all 
sources above the mentioned flux limits at 30 GHz have been 
subtracted from CBI and BIMA data, but the residual signal 
has not been taken into account. We computed the power 
spectrum from residual point sources given the above men- 
tioned flux limits (relative to the two experiments) and the 
two different models. We also took into account the sources 
correlation which contribute about 15-20 per cent. We then 
compared the residual power spectra with the CBI band 
power for I > 2000 and the BIMA ones at I > 5000. Our re- 
sults are presented in fig. El We find that, in the case of the 
emission at the level of model B for all sources, their contri- 
bution would be 2.4% for CBI, and 2.5% for BIMA while, in 
the case of model A the contribution would be 47% to CBI 
and 46% to BIMA. Since Cf z oc al, in the latter case the 
unresolved signal from these sources could produce a bias of 
6 — 7% in the current determination of as with this method 
(Komatsu & Seljak 2002). The real value of as is still a very 
controversial issue, but the one inferred from SZ power spec- 
trum is certainly in the high end of the considered range. A 
lower as value would be in better agreement with determi- 
nations from cluster abundance (Pierpaoli et al. 2003), the 
WMAP results (Spergel et al. 2003) and some weak lensing 
experiments (Heymans et al. 2004). 



3.3 Residual signal in all— sky experiments 

In this section we discuss the residual contribution from 
accretion flow sources to the bispectrum of WMAP and 
Planck. 
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Figure 3. Residual power spectrum at 30 GHz from accretion 
flow sources and the high-1 points from the CBI and BIMA exper- 
iments. Bold-blue (thin-red) lines are relative to the CBI (BIMA) 
flux cut; while the square— blue (triangle— red) points are the CBI 
(BIMA) results with 1 a error bars. The I ~ 9000 BIMA point is 
a 90 per cent upper limit. Dotted (dashed) lines refer to model A 
(B) respectively, as in figure 1. 
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Table 1. Bispectrum values implied by residual accretion-flow 
sources. Column 1 is the frequency of observation, column 2 is the 
assumed flux limit, columns 3 and 4 are bispectrum estimates for 
point source models A and B; column 4 are values reported from 
Argueso et al (2003) for standard radio/IR sources (Toffolatti et 
al. 1998). 



In tableQwe report the bispectrum produced by resid- 
ual sources for different flux cuts and normalizations, com- 
pared with the one implied by more standard radio and in- 
frared (IR) population as in the model by Toffolatti et al. 
(1998) (see also Argueso et al. (2003)) \ WMAP mea- 
sured the bispectrum in Q (40 GHz) and V (61 GHz) 



It should be kept in mind that these quoted values have been 
reduced by a factor 0.8 in order to match the WMAP observations 
(Argueso et al. 2003). 
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bands, and found 6 = (9.5 ± 4.4) x 10~ 5 /iK 3 sr 2 and b = 
(1.1 ± 1.6) x 10~ 5 /^K 3 sr 2 respectively. For a nominal flux 
cut at Sum = 0.75 accretion flow sources would imply 
b = 2.9 x 10 _7 (1.1 x 10~ 8 )AtK 3 sr 2 in the V band for model 
A (B) and b = 1.7 X 10~ 6 (6.0 x l(T 8 )/iK 3 sr 2 in the Q band, 
so they contribute at most at the ***25 per cent level in the 
V band, and 2 per cent in the Q band. Given the WMAP 
bispectrum errors bars, these numbers do not allows us to 
draw constraints on accretion flow sources from these data. 
The current flux limit from WMAP is so high (Su, n — 1 Jy) 
that the contribution from these sources is sub dominant 
with respect to what is inferred by more standard radio/IR 
sources. The observed bispectrum in the current WMAP 
data is consistent with a residual signal produced by the 
population of point sources observed at fluxes S > 1 Jy 
(Komatsu et al. 2003), which is not dominated by accretion 
flow sources, rather by quasars and active galaxies as in the 
Toffolatti et al. (1998) model (see fig-HJ. 

Because of the high sensitivity of Planck, point sources 
are expected to be detected down to a flux between 0.01 and 
0.1 Jy. In tableQwe show the contribution of accretion flow 
sources to the Planck bispectrum. Their contribution at 100 
GHz can be ~ 15 — 20% the one from the Toffolatti et al. 
(1998) radio/IR point sources population as computed in 
Argiieso et al. (2003). In addition, Planck will observe over 
a broad frequency range, allowing a better description of 
point sources properties and possible discrimination between 
accretion-flow sources and other sources. 

3.4 Future high—resolution SZ surveys 

Future high resolution SZ surveys in the 150-350 frequency 
range like ACT or APEX will allow to detect point sources 
down to a flux limit Su m ~ 1 mjy. Because this point source 
population presents an emission that extends to very high 
frequencies (see fig.0, it may give a significant contribution 
to the observed signal in future SZ surveys. Here we report 
the estimated pixel noise produced by residual sources for 
different flux cut and for a fiducial pixel of 1'. Again we 
report our results for the two models shown in Figure 1. 
These two figures should roughly provide a lower and upper 
limit to the contribution coming from elliptical galaxies. Our 
results are summarized in table|2]together with the estimates 
of White & Majumdar (2004) for the population of currently 
observed radio sources. For the sake of comparison, here 
we neglected the contribution of clustering. Accretion flow 
sources may produce a noise that is ~ 25 — 30% of the one 
produced by the population of observed sources. 

At the same frequencies, infrared sources are currently 
predicted to produce a slightly higher signal than radio 
sources (White & Majumdar 2004). These IR estimates, 
however, may be biased high if part of the signal is due 
to lensing effect. Future experiments like SCUBA-2 and the 
Herschel Space Observatory will help in better characteriz- 
ing the IR population. 



4 CONCLUSIONS 

We investigated the possible contribution to current and 
future CMB experiments from the accretion flows around 
supermassive black holes believed to harbor the centers of 



v (GHz) Si im {mJy) a h%A a b _ B cr b>r 



150 1 1.3 0.3 5 (4) 

150 5 2.2 0.5 9 (7) 

150 10 2.8 0.6 12 (9) 

150 50 2.9 0.9 20 (16) 

220 1 1.1 0.3 4 (3) 

220 5 2.0 0.4 7 (5) 

220 10 2.5 0.5 9 (7) 

220 50 2.6 0.7 16 (12) 



Table 2. Pixel noise implied by residual accretion flow sources in 
SZ experiments with 1' resolution. Columns 3 and 4 are predic- 
tions for model A and B. Column 5 are predictions from White 
and Majumdar (2003) on the standard radio sources contribution 
(the two numbers refer to their two different extrapolations from 
low frequencies). The last three columns are in fiK for a nominal 
1' pixel. 



galaxies. These sources are typically faint and numerous, 
therefore they are likely to be observed as a residual signal in 
CMB experiments. A comparison with the number counts of 
WMAP detected sources (and identified as elliptical galax- 
ies) indicates that the bulk of the emission is probably higher 
that the ~ 10 28 erg/s/Hz assumed in our emission model B. 

We have showed that the residual signal produced by 
the emission from these unresolved sources could signifi- 
cantly contribute to the CBI and BIMA observed power 
spectrum. More specifically, these sources could make up 
to ~ 40 — 50% of the observed signal at / > 2000, there- 
fore reducing the inferred value of erf z by 6 — 7%. As for all 
sky experiments, we showed that the residual signal in the 
WMAP maps with the current flux limit is likely to be dom- 
inated by other kinds of sources, with accretion-flow sources 
contributing at the 2-3 per cent level at most. 

Planck will have a better sensitivity than WMAP and 
will allow to detect fainter sources reducing the flux limit 
by at least a factor of 10 and possibly detecting hundreds 
of early-type galaxies. We showed that the residual sig- 
nal from accretion-flow sources may contribute up to the 
15% to the bispectrum around 100 GHz. Unlike other radio 
sources, these sources typically show an inverted spectrum 
and are quite bright also in the infrared, up to about 300 
GHz. This feature may facilitate their detection and dis- 
crimination with respect to the standard radio and infrared 
sources, especially with Planck that offers the potential of an 
all-sky experiment with a broad frequency coverage. A char- 
acterization of the Plank sources with the SDSS will allow 
to clearly distinguish between elliptical galaxies and other 
kinds of sources, allowing to put constraints on accretion 
flow models with radio number counts. 

Future high-resolution CMB experiments like ACT and 
APEX operating in the 150-350 GHz range may also be 
affected by the emission from accretion flows in galaxies. 
The residual power spectrum is predicted to be about 25- 
30% of the one expected for currently observed radio point 
sources. At these frequencies, however, IR sources are likely 
to be dominant, unless their number count is biased high 
due to lensing effects. Future surveys with instruments like 
SCUBA-2 and Herschel will help in clarifying this issue. 
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